Grid re®nement for modeling multiphase ̄ow in discretely fractured porous media
نویسنده
چکیده
A study of the eects of grid discretization on the migration of DNAPL within a discrete-fracture network embedded in a porous rock matrix is presented. It is shown that an insuciently ®ne discretization of the fracture elements can lead to an overprediction of the volume of DNAPL that continues to migrate vertically at the intersection of a vertical and horizontal fracture. Uniform discretization of elements at the scale of one centimetre (or less) accurately resolved the density and capillary pressure components of the head gradient in the DNAPL. An alternative, non-uniform method of discretization of elements within the discrete-fracture network is presented whereby only fracture elements immediately adjacent to fracture intersections are re®ned. To further limit the number of elements employed, the porous matrix elements adjacent to the fracture elements are not similarly re®ned. Results show this alternative method of discretization reduces the numerical error to an acceptable level, while allowing the simulation of ®eldscale DNAPL contamination problems. The results from two ®eld-scale simulations of a DNAPL-contaminated carbonate bedrock site in Ontario, Canada are presented. These simulations compare dierent methods of grid discretization, and highlight the importance of grid re®nement when simulating DNAPL migration problems in fractured porous media. Ó 1999 Published by Elsevier Science Ltd. All rights reserved.
منابع مشابه
Dispersion in Modeling Multiphase Transport through Heterogeneous Geological Media
340 Numerical modeling approaches for multiphase fl ow and tracer or chemical transport in porous media are generally based on methodologies developed for reservoir simulation and groundwater modeling. They involve solving coupled massconservation equations that govern the transport processes of all chemical components in isothermal and nonisothermal subsurface systems using fi nite-difference ...
متن کاملNumerical simulation of non-isothermal multiphase tracer transport in heterogeneous fractured porous media
We have developed a new numerical method for modeling tracer and radionuclide transport through heterogeneous fractured rocks in a non-isothermal multiphase system. The model formulation incorporates a full hydrodynamic dispersion tensor, based on three-dimensional velocity ®elds with a regular or irregular grid in a heterogeneous geological medium. Two dierent weighting schemes are proposed f...
متن کاملA mixed-dimensional finite volume method for multiphase flow in fractured porous media
We present a vertex-centered finite volume method for the fully coupled, fully implicit discretization of two-phase flow in fractured porous media. Fractures are discretely modeled as lower dimensional elements. The method works on unstructured, locally refined grids and on parallel computers with distributed memory. An implicit time discretization is employed and the nonlinear systems of equat...
متن کاملAdaptive Mesh Re nement and Multilevel Iteration for Flow inPorous
An adaptive local mesh re nement algorithm originally developed for unsteady gas dynamics by Berger [12, 14] is extended to incompressible ow in porous media. Multilevel iteration and domain decomposition methods are introduced to accommodate the elliptic/parabolic aspects of the ow equations. The algorithm is applied to a two-phase polymer ooding model consisting of a system of nonlinear hyper...
متن کاملSmoothed particle hydrodynamics and its applications for multiphase flow and reactive transport in porous media
Smoothed particle hydrodynamics (SPH) is a Lagrangian method based on a meshless discretization of partial differential equations. In this review, we present SPH discretization of the Navier-Stokes and advection-diffusionreaction equations, implementation of various boundary conditions, and time integration of the SPH equations, and we discuss applications of the SPH method for modeling pore-sc...
متن کامل